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Summary 

Two Euler analysis techniques— finite difference and finite 
volume— were employed to predict the blade surface pressure 
distributions of a large-scale advanced propeller. The predicted 
pressure distributions were compared with wind tunnel data. 
Both techniques produced blade pressure distributions that are 
in fairly good agreement with the data over the range of test 
Mach numbers of 0.2 to 0.78. However, the numerical 
simulations fail to predict correctly the measured pressure 
distributions for the low Mach number, high-power case. The 
data indicate the presence of a leading-edge vortex for this 
case. A discussion of the compressibility effects is also 
presented. 

Introduction 

Advanced propellers employ highly loaded, highly swept, 
thin blades with variable pitch to achieve a higher propulsive 
efficiency than is achieved by the current high-bypass turbo- 
fans. Model tests have been carried out in wind tunnels and 
in flight to understand the aerodynamics and acoustics of the 
advanced designs. The propeller blade flows are complicated 
and involve a wide range of flow regimes, such as (1) takeoff 
conditions, where leading edge and tip vortices influence the 
performance and noise generation; and (2) cruise conditions, 
where shock waves add to the complexity. An understanding 


of the complex advanced propeller blade flows is essential for 
improving the current design methodologies. 

Scale-model tests have provided integrated performance 
data. To obtain aerodynamic blade data, a series of tests on 
a large-scale 9-ft single-rotation SR7L propeller were con- 
ducted in a transonic wind tunnel in Modane, France. The 
propfan was operated in a two-bladed configuration (fig. 1). 
This configuration was adopted for the tests because the drive 
power was insufficient for the eight-bladed propeller. The 
design characteristics of the SR7L propeller are given in 
table I. 

The blade surface pressure was measured on a specially 
designed, pressure-tapped SR7L blade. There were 13 radial 
stations of pressure taps on each surface of the blade: 16 
pressure taps were distributed along the blade chord at each 
radial station on the pressure (face) side and 20 were distributed 
on the suction (camber) surface. The taps were concentrated 
on the blade surface in the areas where large pressure gradients 
were expected. Each chordwise distribution was aligned 
approximately with the expected streamline on the blade at 
the design cruise condition (fig. 2). The blade surface pressures 
were recorded one radial station at a time, and the tests were 
conducted over the Mach number range of 0.03 to 0.78. A 
description of the nature and extent of the large-scale advanced 
propeller wind tunnel tests is given by Campbell et al. (ref. 1). 

The purpose of the present investigation is to evaluate the 
ability of the Euler analysis techniques to predict the blade 



surface pressure distributions. The three-dimensional Euler 
equations are solved by two different analysis techniques: (1) 
a finite difference technique (NASPROP) originally developed 
by Bober et al. (ref. 2), and (2) a finite volume technique 
developed by Denton (ref. 3). In this paper, the predicted 
chordwise pressure distributions and the integrated power 
coefficients are compared with the test data for the Mach 
number range from 0.2 to 0.78. 

Computational Methods 

Two numerical codes are employed to simulate the flow field 
about the two-bladed single-rotation propeller (SR7L) subject 
to conditions prescribed in the Modane, France, test. Both 
codes use an asymptotic time-marching approach to obtain the 
steady-state solutions of the three-dimensional Euler equations 
governing the propeller flow. Brief descriptions of each code 
follow. 

Finite Difference Euler Solver (NASPROP) 

The finite difference Euler solver (NASPROP) was 
originally developed by Bober et al. (ref. 2). Their three- 
dimensional Euler analysis was the first to provide detailed 
solutions of the advanced propeller flow field. It employs the 
numerical method of Beam and Warming. Since then the code 
has been rewritten incorporating improved boundary condi- 
tions and has been vectorized for the CRAY computer. 

The three-dimensional, rotational, inviscid flow field about 
a propeller is governed by the Euler equations. The equations, 
in conservative differential form, are transformed from a 
cylindrical coordinate system to a time-dependent, body-fitted, 
curvilinear reference frame. The transformed equations are 
treated as a system of coupled first-order, autonomous differen- 
tial equations. The equations are cast in nondimensional form. 
First-order-accurate implicit Euler time differencing is 
employed, and the flux vectors are linearized by Taylor series 
expansion up to second order in time. The spatial derivatives 
are replaced by central differences with appropriate one-sided 
differences at the boundaries. The resulting implicit operators 
involve inversions of a block tridiagonal system. The system 
can be reduced to scalar tridiagonal matrices by applying a 
diagonalization technique (ref. 4) to the flux Jacobian matrices. 
Constant coefficient, second-order implicit and fourth-order 
explicit dissipation terms are added for stability considerations. 
A periodic tridiagonal matrix inversion is employed in the 
circumferential direction. Further details of the solution 
procedure and the implementation of the boundary conditions 
can be found in references 5 and 6. 

Finite Volume Euler Solver (Denton) 

The finite volume Euler solver was originally developed by 
Denton for the analysis of turbomachinery flows (ref. 3) and 
was later modified to compute the advanced propeller flow 


field (refs. 7 and 8). In Denton’s analysis, the three- 
dimensional Euler equations are written in a coordinate frame 
rotating with the blade. The governing equations are solved 
in the physical domain on a grid that is generated algebraically. 
The flow domain is discretized by using three kinds of sur- 
faces: the bladewise surface, the streamwise surface, and the 
spanwise surface. The computational nodes are located at the 
comers of each volume element. The fluxes of mass, momenta, 
and energy through each face are calculated by using the 
averages of flow properties stored at the comers of that face. 

Denton’s explicit time-marching method employs an opposed 
difference scheme. The scheme uses upwind differences for 
fluxes of mass and momenta, and downwind differences for 
pressure in the streamwise direction. The derivatives in the 
circumferential direction are evaluated by using central 
differences. The time-marching method uses a time step such 
that the Courant number is close to unity. The method itself 
ensures stability by solving the flow equations in the order 
of continuity, energy, and momentum. A variable time step 
appropriate for each volume element is used to accelerate 
convergence. The method is of first-order accuracy and is 
modified to achieve greater accuracy by adding a lagged 
correction factor to correct the down winded pressure. A 
smoothing factor is used after each time step to smooth out 
any waviness in the circumferential direction. The convergence 
is judged on the basis of mass conservation and the maximum 
change in the axial velocity component. Further details of the 
solution procedure and the implementation of the boundary 
conditions can be found in references 3 and 8. 


Computational Grid 

In this paper, the steady-state propeller flow field with the 
propeller axis at zero angle of attack is computed. Hence only 
one blade passage needs to be considered as a computational 
domain since the flow field is periodic in the blade-to-blade 
direction. 

The computational domain is discretized by three types of 
surfaces. The bladewise surfaces are bounded by two planes— 
one plane containing the suction side of one blade and the other 
plane containing the pressure side of the adjacent blade 
(fig. 3(a)). The surfaces are clustered near the blades, and 
larger spacings are used between the blades. The streamwise 
surfaces (fig. 3(b)) are placed in the radial direction and are 
limited by the spinner-nacelle surface and the far field 
boundary. The far field boundary is placed at two blade radii 
beyond the blade tip. Small spacings are used in the blade tip 
region to capture the tip vortex flow. The spanwise surfaces 
are bounded by inlet and exit planes in the axial direction 
(fig. 3(c)). The inlet plane is placed at one radius upstream 
of the spinner, and the exit plane is located at three radii behind 
the blade. Grid points are clustered near the leading and trailing 
edges. 
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From a few grid-sensitivity calculations, 35 chord wise and 
20 spanwise points on each side of the blade are found to be 
adequate to obtain a solution that is nearly independent of the 
grid. Figure 3(d) shows the distribution of grid points on the 
blade surface. A total of 61 by 33 by 35 (axial, radial, and 
circumferential) grid points is used to discretize the whole 
computational space. The same grid density on the blade 
surfaces, and similar grid arrangements, are employed by the 
two computer programs. 


Results and Discussion 

The steady three-dimensional flow field of an advanced 
propeller, SR7L, is considered here. The Euler equations 
governing the flow field are solved by a finite difference and 
a finite volume technique. The numerical results are compared 
with wind tunnel data. 

The numerical data for nine of the Modane test cases 
(table II) were generated by the two numerical codes 
independently. The data input to the codes were the same. In 
particular, exactly the same blade setting angles were used in 
both codes. Although each code used its own grid generator 
for high Mach number runs, both used the same grid generator 
for low Mach number (0.2) runs. The results are compared 
with the wind tunnel data in the form of surface pressure 
contours, chordwise pressure distributions, and total power 
coefficients. 

The pressure coefficient is defined in this paper as 


Cp = V 

U 


Pjj ~ Poo 

0-5poo(F^ + rjjw 2 ) 


where Py and r t j are the local pressure and the radial distance 
to the data point; p*,, and V » are the pressure, density, 
and velocity of the free stream; and co is the rotational speed 
(rad/sec) of the blade. The computed pressure is interpolated 
to match the radial location of the pressure taps. The total 
power coefficient of the propeller is computed from the blade 
surface pressure by 


CP E P Jl 2 D i r ^ P U A ^ij)eB 
ij 

where A A y is the elemental surface area formed by four 
neighboring grid points, r {j is the radial distance of the 
element, D is the diameter of the propeller, fi is the rotational 
velocity (rev/sec), and B is the number of blades. The subscript 
0 denotes the component in the tangential direction. 

The Modane test results (ref. 9) and the computed power 
coefficients are presented in table n. The Mach number and 
the advance ratio for the computations were matched with the 
experimental conditions. The use of nominal blade setting 


angles for computations resulted in power coefficients higher 
than the measured ones. Based on the computed power 
coefficients, a second set of blade setting angles were chosen 
within the tolerance limits, to make the predictions closer to 
the data. These blade setting angles, which were used in the 
computations, and the resulting power coefficients are given 
in the table. For high Mach number tests, the power 
coefficients are still overpredicted. The blade setting angles 
may have to be lowered beyond the tolerance limits to improve 
the agreement in these cases. No further change in the blade 
setting angle was made. 

The power coefficients predicted by the numerical codes 
closely follow the power trends that are observed in the 
experiments. Tests 5 and 6 are designed to compare the power 
at cutback and takeoff conditions at low Mach number. The 
predictions show the correct trend, except that the NASPROP 
code underpredicts the measured value in test 6. 

Tests 7, 8, and 9 show the power trend when the blade angle 
is changed while the Mach number and the advance ratio are 
kept constant. For low-power range, NASPROP results show 
excellent agreement with the wind tunnel data. 

Tests 10 through 12 are designed to show the Mach num- 
ber effect when the blade angle is kept constant. Through 
these tests the power stays relatively constant, although the 
results of the finite difference code show a slight increase with 
Mach number. The Denton code results show nearly the same 
trend as the data. The last test was conducted as a variation 
of test 12 with a different advance ratio. The inverse relation 
between the power coefficient and the advance ratio is correctly 
predicted by both numerical codes; however, the finite 
difference code tends to overpredict the power consistently 
for these cases. 

Three representative test cases were chosen for detailed 
comparisons: test 12 for a high Mach number (0.775), test 8 
for an intermediate Mach number (0.5), and test 5 for a low 
Mach number (takeoff, 0.2). 

Contours of measured surface pressure for test 12 are shown 
in figure 4. Because of the coarseness of the experimental data 
points, the blade shape and the contour lines are not smoothly 
constructed. Although there are 20 pressure taps in each of 
the 13 radial stations on the suction side and 16 taps on each 
pressure surface radial station, data were recorded for only 
10 radial stations for test 12. Also, a few pressures were not 
valid at several radial stations. Missing data were approximated 
by interpolation or extrapolation from neighboring data. 

The corresponding data generated by the two numerical 
codes are also shown in figure 4. Similarity between the 
numerical results is clearly seen. On the suction side (fig. 4(a)), 
the contour features can be easily identified with those of the 
data. The agreement of the predictions with the data on the 
pressure surface (fig. 4(b)) is not as good as that on the suction 
surface. 

Detailed chordwise pressure distributions at four radial 
stations (i r/R = 0.569, 0.739, 0.862, and 0.975) are shown 
in figure 5(a), (b), (c), and (d), respectively. Better agreement 


3 



is seen on the suction side, except that a strong shock is 
generated at the outboard trailing edge by the NASPROP code. 

Although some test data are missing in this region and it 
is difficult to assess exact flow, a shock as strong as that 
predicted by NASPROP may not have developed during the 
test. Contrary to the NASPROP results, the Denton code does 
not show the formation of a strong shock though it clearly 
shows the compressibility effect in the region and agrees well 
with the measurements. Both codes tend to overpredict the 
pressures on the pressure side, which is reflected in the power 
coefficient comparison in table n. The inverted loading near 
the leading-edge region obtained in the data is predicted poorly 
on the inboard region. This may be the result of the leading- 
edge geometry represented by the numerical codes not exactly 
matching the physical blade geometry. The estimated structural 
deflections at the design point cruise condition were used for 
all of the runs. 

The surface pressures produced by the NASPROP code are 
compared with the test 8 data in figure 6, and the local 
chordwise pressure distributions at the four radial stations are 
compared with test 8 data in figure 7. As in the previous case, 
good agreement is shown on the suction side. Except on the 
blade edges, the contour levels on the pressure side compare 
reasonably well. It should be noted that the blade outline in 
the test data is based on the pressure tap locations and does 
not represent the true outline of the blade. In general, better 
agreements are seen in the outboard section of the blade, and 
discrepancies are observed at the inboard leading-edge section. 

Chordwise pressure profiles (fig. 7) reveal that the loading 
is fairly evenly distributed from the hub to the tip at this 
condition. Both codes show reasonably good agreement at all 
four radial stations. The pressure side is predicted better by 
the Denton code than by the NASPROP code. 

Numerical results are compared with experimental test 5 data 
(Mach 0.2) in figures 8 and 9. The agreement of the predicted 
pressure contours with data (fig. 8) is good except for a few 
isolated places like the edges of the blade and the hub region. 
Figure 9 shows that both numerical predictions agree very well 
with test data, especially in the outboard section of the blade. 
The Denton code slightly overpredicted the pressure loading 
throughout the blade, which resulted in a power coefficient 
that was higher than the experimental value (table II). 

Figure 10 shows a low Mach number (0.2) case with a higher 
power (test 6), and presents an interesting flow phenomenon. 
The numerical results are represented by the Denton 
calculations. The code is able to capture the high-pressure 
loading at the leading edge reasonably well. The wind tunnel 
data show broadening of the suction peak on the suction side 
at radial stations r/R = 0.862 and r/R = 0.905. This behavior 
of the pressure curve is considered to be the result of the 
formation of a leading-edge vortex. Oil flow studies of surface 
streamlines indicate the formation of such leading-edge 
vortices (refs. 10 and 1 1). Unfortunately, with the conditions 
prescribed in the Modane test and the grid distribution 
employed in the present numerical solutions, the vortex flow 


has not been captured by either numerical code. Grid refine- 
ment and improved dissipation models may be necessary to 
predict the leading-edge vortex. Numerical solutions of the 
Euler equations have been found to predict the leading-edge 
vortices on swept wings fairly well (ref. 12). 

Lastly, the Mach number effect is demonstrated in figures 1 1 
and 12. Figure 1 1(a) shows the NASPROP pressure contours 
on the suction side for four Mach numbers, 0.5, 0.6, 0.7, and 
0.775 (tests 8, 10, 11, and 12, respectively) and the corre- 
sponding experimental data. The pressure surface predictions 
are compared with data in figure 1 1(b). It should be noted that 
for this compressibility study the blade angle is kept constant. 
However, for the experimental results, a slightly different 
blade angle was used for the Mach 0.5 case (see table II). With 
increasing Mach number, more significant changes in the flow 
occur on the suction side than on the pressure side. The 
predicted flow features compare well with the experimental 
data. 

Figure 12 shows the local pressure distribution at two out- 
board stations, r/R — 0.806 and 0.939, for the compressibility 
study. Although no significant changes are observed for low 
Mach numbers, the evidence of a compression wave starts to 
develop on the suction side when the Mach number reaches 
about 0.7. At a free-stream Mach number of 0.78, the com- 
pression wave fully develops to a trailing-edge shock. The 
NASPROP results agree with the experimental data up to the 
development of the compression wave. However, the strong 
trailing-edge shock formation is not indicated by the test data 
or the results of the Denton code. The reasons for this 
discrepancy are not clear. 


Concluding Remarks 

The detailed chordwise pressure distributions of an advanced 
propeller were predicted by two numerical techniques of 
solving the three-dimensional Euler equations— finite 
difference and finite volume— and were compared with wind 
tunnel data. In general, both techniques predicted the blade 
surface pressure distributions fairly well for the range of Mach 
numbers considered (0.2 to 0.78). The strong trailing-edge 
shock predicted by the finite difference code (NASPROP) at 
high Mach numbers (0.775 and 0.78) was not observed either 
in the experimental data or in the Denton code predictions. 
In the intermediate Mach number range (0.5 to 0.7) the 
agreement of the predictions with data was good. Generally, 
the predictions showed better agreement at the outboard section 
than at the inboard section. For the low Mach number (0.2) 
cases, the agreement of the predictions with data was best for 
the low-power case. The data for the high-power case seems 
to show the formation of a leading-edge vortex that is not 
captured by the numerical techniques. Grid refinement and 
improved dissipation models may be required to capture the 
leading-edge vortex. 
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TABLE I.-SR7L PROPELLER 
DESIGN PARAMETERS 


Disk loading, kW/m 2 (shp/diam 2 ) 259.6 (32) 

Tip speed, m/sec (ft/sec) 243.8 (800) 

Tip diameter, m (ft) 2.743 (9) 

Hub-to-tip ratio 0.24 

Advance ratio 3.06 

Power coefficient 1.448 

Integrated design lift coefficient 0.191 

Blade tip sweep 36 

Number of blades 8 


TABLE n.-MODANE TEST OPERATING CONDITIONS 


Test 

number 

Mach 

number 

Advance 

ratio 

Blade angle, 
deg 

Power coefficient 

Test 

NASPROP 

calculation 

Denton 

calculation 

Test 

Calculation 

5 

0.2 

0.88 

26. 6± 1.0 

25.6 

0.100±0.002 

0.108 

0.121 

6 

.2 

.883 

31.3 ± .9 

30.4 

.250 ±.001 

.193 

.283 

7 

.5 

3.083 

58.5 ± 1.0 

57.5 

.642 ±.009 

.564 

.692 

8 

.5 

3.067 

55.1 ± .8 

55.0 

.361 ±.005 

.369 

.344 

9 

.5 

3.071 

51.8 ± .9 

51.6 

.110±.004 

.106 

.112 

10 

.6 

3.078 

54.5 ± 1.3 

53.6 

.230 ±.006 

.254 

.262 

11 

.7 

3.064 

54.5 ± 1.3 

53.6 

.228 ±.006 

.271 

.244 

12 

.775 

3.088 

54.5 ± 1.3 

53.6 

.227 ±.003 

.286 

.247 

13 

.78 

3.209 

54.5 ± 1.3 

53.6 

. 1 1 1 ± .005 

.207 

.134 
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Figure 1.— SR7L propfan installed in the Sl-MA transonic test section. 



Figure 2.— Pressure tap locations. 
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(d) 


(a) Blade wise surfaces. 

(b) Stream wise surfaces. 

(c) Span wise surfaces. 

(d) Distribution of grid points on the blade surface. 
Figure 3.— Typical computational grid. 
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(a) Suction surface. 

(b) Pressure surface. 

Figure 4.— Comparison of predicted and measured pressure coefficient distributions for test 12 (Mach 0.775). 
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FRACTION OF CHORD, X/C 

(a) Radial location, r/R , 0.569. 

(b) Radial location, r/R , 0.739. 

(c) Radial location, r/R , 0.862. 

(d) Radial location, r/R , 0.975. 

Figure 5.— Chordwise distribution of pressure coefficients for test 12 (Mach 0.775). 
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(a) 



NASPROP PREDICTION 


(b) C-8S-061W 

(a) Suction surface. 

(b) Pressure surface. 

Figure 6.— Comparison of predicted and measured pressure coefficient distributions for test 8 (Mach 0.5). 
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PRESSURE COEFFICIENT, Cp 


EXPERIMENTAL DATA 
NASPROP PREDICTIONS 
DENTON PREDICTIONS 





(a) Radial location, r/R , 0.569. 

(b) Radial location, r/R , 0.739. 

(c) Radial location, r/R , 0.862. 

(d) Radial location, r/R , 0.975. 

Figure 7.— Chordwise distribution of pressure coefficients for test 8 (Mach 5). 
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(b) 0-83-061*15 

(a) Suction surface. 

(b) Pressure surface. 


Figure 8.— Comparison of predicted and measured pressure coefficient distributions for test 5 (Mach 0.2). 
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(a) Radial location, r/R, 0.569. 

(b) Radial location, r/R , 0.739. 

(c) Radial location, r/R , 0.862. 

(d) Radial location, r/R , 0.975. 

Figure 9.— Chordwise distribution of pressure coefficients for test 5 (Mach 0.2). 
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Figure 11.— Comparison of predicted and measured compressibility effects: pressure coefficient distributions for various Mach numbers. 
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(b) Pressure surface. 
Figure 1 1.— Concluded. 
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Figure 12.— Comparison of predicted and measured compressibility effects: chordwise distribution of pressure coefficients for two radial locations for various 
Mach numbers. 
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